TIME REVERSAL OF VOLTERRA PROCESSES DRIVEN 
STOCHASTIC DIFFERENTIAL EQUATIONS 



L. DECREUSEFOND 



Abstract. We consider stochastic differential equations driven by some Volterra 
processes. Under time reversal, these equations are transformed into past de- 
pendent stochastic differential equations driven by a standard Brownian mo- 
tion. We are then in position to derive existence and uniqueness of solutions 
of the Volterra driven SDE considered at the beginning. 



1. Introduction 

Fractional Brownian motion is one the first example of a process which is not a 
semi-martingale and for which we aim to develop a stochastic calculus. That means 
we want to define a stochastic integral and solve stochastic differential equations 
driven by such a process. From the very beginning of this program, two approaches 
do exist. One approach is based on the sample-paths properties of fBm, mainly 
its Holder continuity or its finite p-variation. The other way to proceed relies on 
the gaussiannity of fBm. The former is mainly deterministic and was initiated 
by Zahle [41], Feyel, de la Pradelle [12] and Russo, Vallois [31, 32]. Then, came 
the notion of rough paths introduced by Lyons [22], whose application to fBm 
relies on the work of Coutin, Qian[4]. These works have been extended in the 
subsequent works [20, 21, 3, 15, 14, 23, 26, 27, 16, 25, 8]. A new way of thinking 
came with the independent but related works of Feyel, de la Pradelle [13] and 
Gubinelli [17]. The integral with respect to fBm was shown to exist as the unique 
process satisfying some characterization (analytic in the case of [13], algebraic in 
[17]). As a byproduct, this showed that almost all the existing integrals throughout 
the literature are all the same as they all satisfy these two conditions. Behind each 
approach but the last too, is a construction of an integral defined for a regularization 
of fBm, then the whole work is to show that under some convenient hypothesis, the 
approximate integrals converge to a quantity which is called the stochastic integral 
with respect to fBm. The main tool to prove the convergence is either integration by 
parts in the sense of fractional deterministic calculus, either enrichment of the fBm 
by some iterated integrals proved to exist independently or by analytic continuation 
[37, 36]. 

In the probabilistic approach [7, 6, 5, 9, 19, 30, 29, 2, 1], the idea is also to define 
an approximate integral and then prove its convergence. It turns out that the key 
tool is here the integration by parts in the sense of Malliavin calculus. 

In dimension greater than one, with the deterministic approach, one knows how 
to define the stochastic integral and prove existence and uniqueness of fBm driven 
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SDEs for fBm with Hurst index greater than 1/4. Within the probabilistic frame- 
work, one knows how to define a stochastic integral for any value of H but one 
can't prove existence and uniqueness of SDEs whatever the value of H. The pri- 
mary motivation of this work was to circumvent this problem. 

In [7, 9], we defined stochastic integrals with respect to fBm as a "damped- 
Stratonovitch" integral with respect to the underlying standard Brownian motion. 
This integral is defined as the limit of Ricmann-Stratonovitch sums, the convergence 
of which is proved after an integration by parts in the sense of Malliavin calculus. 
Unfortunately, this manipulation generates non-adaptiveness: Formally the result 
can be expressed as 

t 

u(s) o dB H (s) = 6(K,* t u) + trace(/C*Vu). 

Even if u is adapted (with respect to the Brownian filtration), the process (s i-> 
K.fU(s)) is anticipative. However, the stochastic integral process (t i-> j u(s) o 
dB H (s)) remains adapted so the anticipativeness is in some sense artificial. The 
motivation of this work is to show that up to time reversal, we can work with 
adapted process and Ito integrals. In what follows, there is no restriction about the 
dimension but we need to assume that for any component B H is an fBm, the Hurst 
index of which is greater than 1/2. 

Consider that we want to solve the equation 

(1) X t =x+ f a(X s ) o dB H (s), < t < T 

Jo 

where a is a deterministic function whose properties will be fixed below. It turns 
out that it is essential to investigate the more general equations: 

(A) X r>t = x+ ( a(X TiS ) o dB H (s), < r < t < T. 

J r 

The strategy is then the following: We will first consider the reciprocal problem: 

(B) Y r>t = x- [ <r(Y Sit ) o dB H (s), < r < t < T. 

J r 

The first critical point is that when we consider {Z r ^ t ■— Y t - r ^ t , r G [0, t]}, this 
process solves an adapted, past dependent, stochastic differential equation with 
respect to a standard Brownian motion. Moreover, because Kh is lower-triangular 
and sufficiently regular, the trace term vanishes in the equation defining Z. We 
have then reduced the problem to an SDE with coefficients dependent on the past, 
a problem which can be handled by the usual contraction methods. This paper is 
organized as follows: After the preliminaries of Section 2, we address, in Section 
3, the problem of Malliavin calculus and time reversal. This part is interesting 
in its own since stochastic calculus of variations is a framework oblivious to time. 
Constructing such a notion of time is achieved using the notion of resolution of 
the identity as introduced in [40]. We then introduce the second key ingredient 
which is the notion of strict causality or quasinilpotence, see [42] for a related 
application. In Section 4, we show that solving Equation (B) reduces to solve a 
past dependent stochastic differential equation with respect to a standard Brownian 
motion, see Equation (C) below. In Section 5, we prove existence, uniqueness and 
some properties of this equation. Technical lemmas are postponed to Section 6. 
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2. Preliminaries 

Let T > be fix real number. For a function / G £ 1 ([0, T]; R"), we define r T f 

by 

t t/(s) = f(T - s) for any S £[0, T]. 
For t G [0, T], e t f will represent the restriction of / to [0, t], i.e., e t f = /lro,ti- For 
any linear map A, we denote by A^, its adjoint in £ 2 ([0, T]; R"). For 77 G (0, 1], 
the space of 77-Holder continuous functions on [0, T] is equipped with the norm 

II/IIhoi ( „)= sup m * ) ~ f , i " )l + \\f\U 

0<s<t<T \t - S\ I 

Its topological dual is denoted by Hol(r?)*. For / G £*([(), T]; R n ; dt), (denoted by 
C 1 for short) the left and right fractional integrals of / are defined by : 



= YmJ /(*)(* - dt , x <T, 

where 7 > and 7°+ = 7°_ = Id . For any 7 > 0, p, q > 1, any / G £ p and g E C q 
where + < 7, we have : 

(2) f f(s)(r 0+ g)(s) ds = f(P T -f)( S )g{ S ) ds. 
Jo Jo 

The Besov-Liouville space ij+ (£ p ) := Z+ p is usually equipped with the norm : 

(3) K 7 +/llz+ p = ll/lk- 



Analogously, the Besov-Liouville space (£ p ) := I is usually equipped with 
the norm : 

ll^/llz-, = ll/IU- 
We then have the following continuity results (see [12, 33]) : 

Proposition 2.1. i. If < 7 < 1, 1 < p < I/7, then 7j + is a bounded operator 

from C p into C q with q = p(l — ~/p)~ 1 . 
ii. For any < 7 < 1 and any p > 1, Z+ p is continuously embedded in Hol(7 — 1/p) 

provided that 7 — l/p > 0. 
Hi. For any < 7 < ft < 1, Hol(/3) is compactly embedded in I 1 >oa . 
iv. For "fp < 1, the spaces Z+ p and X~ p are canonically isomorphic. We will thus 

use the notation X 7;P to denote any of this spaces. 



3. Malliavin calculus and time reversal 

Our reference probability space is Q = Co([0,T], R n ), the space of R"-valued, 
continuous functions, null at time 0. The Cameron-Martin space is denoted by 
H and is defined as H = l£ + (£ 2 ([0, T])). In what follows, the space £ 2 ([0, T]) is 
identified with its topological dual. We denote by k the canonical embedding from 
H into fl. The probability measure P on £1 is such that the canonical map W : uj 
(u>(t), t G [0, T]) defines a standard n-dimensional Brownian motion. A mapping 
(f> from £1 into some separable Hilbert space ^ is called cylindrical if it is of the 
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form 4>{w) = J2i=i M.i v is,w), • • • , (vi >n ,w))xi where for each i, E C~(R™,R) 
and (vij , j = 1 . . . n) is a sequence of ft*. For such a function we define V w </> as 

X/ W (t>{w) = ^ 9jfi{(Vi,l, w), ■ ■ ■ ,{Vi, n ,w))Vi,j ®Xi, 

where v is the image of v G ft* by the map (7q + ok)*. From the quasi- invariance 
of the Wiener measure [39], it follows that V w is a closable operator on L p (£l;f)), 
p > 1, and we will denote its closure with the same notation. The powers of V w 
are defined by iterating this procedure. For p > 1, fc £ N, we denote by IDp ;fc (.f)) 
the completion of i>valucd cylindrical functions under the following norm 



ii0iip, fe = Eii( vW )>n^^ 



®£p([0,1])®*) ■ 



i=0 



We denote by L p4 the space D Pj i(/7([0, T]; R"). The divergence, denoted <5 W is 
the adjoint of V w : v belongs to Dom p (5 W whenever for any cylindrical </>, 



and for such a process v, 



E 



/ Vs 

Jo 

/ Vs^ 

Jo 



ds 



< c 



I Lp 



ds 



= E\6S w v] 



We need first to introduce the "time reversal" operator, denoted by tt and defined 
by: 



r T : £°([0,T];R") 



W (T-.)- 



We introduced the temporary notation W for standard Brownian to clarify the 
forthcoming distinction between a standard Brownian motion and its time reversal. 
Actually, the time reversal of a standard Brownian is also a standard Brownian 
motion and thus, both of them "live" in the same Wiener space. We now precise 
how their respective Malliavin gradient and divergence are linked. Consider B = 
(B(t), t € [0, T]) an n-dimcnsional standard Brownian motion and B T = (B(T) — 
B(T — t),t& [0, T]) its time reversal. Consider the following map 



6 T : ft 



to - 



bj(T) — T T U!, 



and the commutative diagram 



0+ . 
ft D H 



Heft 
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Note that 6^ = 6 T 


since w(0) = 


: 0. For a function / G C fc °°(R™ fc ), we define 


V r /(w(ti),- 


• ' , w(tfe)) = 


k 

I>/Mfi)>-- 


• , w(ife))l[o,t 3 ]W and 


V r /(di(tl),- 


• ' , w(ife)) = 


Ew(ti).-- 





The operator V = V s (respectively V = V s ) is the Malliavin gradient associated 
with a standard Brownian motion (respectively its time reversal). Since, 

/(£>(*i), • • • , oj(t k )) = f(w(T) - lu(T - h), • • • , w(T) - oj(T - t k )), 

we can consider f(u)(ti), • • • , £>(ifc)) as a cylindrical function with respect to the 
standard Brownian motion. As such its gradient is given by 

k 

V r /(w(ti), • • • , w(tfe)) = • • • ,w(^))l [T _ tjiT] (r). 

We thus have, for any cylindrical function F, 

(4) VFoGtH =t t VF(w). 

Since O^P = P and t t is continuous from C p into itself for any p, it is then easily 
shown that the spaces D Pj k and D P: & (with obvious notations) coincide for any p, k 
and that (4) holds for any element of one of theses spaces. Hence we have proved 
the following theorem: 

Theorem 3.1. For any p > 1 and any integer k, the spaces D Pj k and Pj k coincide. 
For any F e D P) k f or some p, k, 

V(F o 6 T ) = t t V(F o 9 T ), P a.s.. 

By duality, an analog result about follows for divergences. 

Theorem 3.2. A process u belongs to the domain of 8 if and only if t t u belongs 
to the domain of 5 and then, the following equality holds: 

(5) 5(u(lu))(lu) = S(t t u(oj))(lu) = 8{t t u o 9 r )(w). 
Proof. For fee£ 2 , for cylindrical F, we have on the one hand: 



E 

and on the other hand, 
E 



F(uj)5h(u>) = E (VF(cD), h) C 2 



(VF(w), h)c\ = E 1(t t VF o 6tM, h)cA 
= E[(VFo9 T (w), t t /i) £2 ] 

= E[foe T (^(T T ft)H] 

= E[F(w)5(T T h)(oj)} . 

Since this is valid for any cylindrical F, (5) holds for h € C 2 . Now, for tt in the 
domain of divergence (see [28, 39]), 

5u = y^((u, h^&Shi - (Vu, ft» ® hi) C 2^ C 2^ , 



6 



Time reversal of fBm driven SDEs 



where (hi, i S N) is an orthonormal basis of £ 2 ([0, T]; R"). Thus, we have 
S(u(lu))(uj) = y^ y ((u(u}), hi) C 2Shi(u}) - (Vu(w), hi ® /i;)/^^^ 

= ^((u(o>), hi) C 25(T T hi){uj) - (Vu(u),T T h t <8> hi) C 2® C 2^j 

i 

= y^f(T r M(a;), T T hi) C 2S(T T hi)(uj) - (Vr T u(w), r T /ij <g> T T hi)c^^c 2 ^ , 

where we have taken into account that r T in an involution. Since (/ij, i e N) is an 
orthonormal basis of £ 2 ([0, T]; R n ), identity (5) is satisfied for any u in the domain 
of 5. □ 

3.1. Causality and quasi-nilpotence. In anticipative calculus, the notion of 
trace of an operator plays a crucial role, we refer to [10] for more details on trace. 

Definition 3.1. Let V be a bounded map from £ 2 ([0, T]; R") into itself. The map 
V is said to be trace-class whenever for one CONB (h n , n > 1) of £ 2 ([0, T}; R"), 

^2 \ (Vh n , h n ) C 2\ is finite. 

n>l 

Then, the trace of V is defined by 

tracc(V) = y^(Vfo n , h n ) C 2. 

n>l 

It is easily shown that the notion of trace does not depend on the choice of the 
CONB. 

Definition 3.2. A family E of projections E\, A e [0, 1] in a Hilbert space H is 
called a resolution of the identity if it satisfies the conditions 

(1) E Q = and E 1 = Id. 

(2) E\E^ = E\ Afl . 

(3) lim^A Ef, = E x for any A G [0, 1) and lim Mt i = Id . 

For instance, the family £J = (cat, A e [0, 1]) is a resolution of the identity in 
£ 2 ([0, T]; R"). 

Definition 3.3. A partition ir of [0, T] is a sequence {0 = to < t\ < . . . < t n = T}. 

Its mesh is denoted by \ir\ and defined by \ir\ — swp i \ti+i — U\. For t G 7r\{T} 7 t+ 
is the least term of it strictly greater than t. 

The causality plays a crucial role in what follows. The next definition is just the 
formalization in terms of operator of the intuitive notion of causality. 

Definition 3.4. A continuous map V from an Hilbert space H into itself is said 
to be E-causal if the following condition holds: 

E X VE X = ExV for any A £ [0, 1]. 

For instance, for H = £ 2 ([0, T); R"), an operator V in integral form Vf(t) = 
J Q T V(t, s)f(s) ds is causal if and only if V(t, s) = for s > t, i.e., computing V f(t) 
needs only the knowledge of / up to time t and not after. Unfortunately, this notion 
of causality is insufficient for our purpose and we are led to introduce the notion of 
strict causality as in [11]. 
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Definition 3.5. Let V be a causal operator. It is a strictly causal operator whenever 
for any e > there exists a partition ir of [0, T] such that for any it' C it, 

\\(E t , + - E t ,)V{E t , + - E t ,)\\ < e for any t' e rr'. 

Note carefully that the identity map is causal but not strictly causal. Indeed, if 
V = M, 

\\{E t , + - E t ,)V(E t , + - E t ,)\\ = \\E t , + - E t ,\\ = 1 
since E t i — E v is a projection. However, for 7 > 0, we have the following result: 

Lemma 3.3. Let H = £ 2 ([0, T); R"). Assume the resolution of the identity to be 
either E = (e\T, A £ [0, 1]) or E = (Id-e( 1 _ A ) T , A e [0, 1]). IfV is an E-causal 
map continuous from C 2 into C p for some p > 2 then V is strictly E-causal. 

Proof. Let tt be any partition of [0, T]. Assume E = (e\T, A e [0, 1]), the very 
same proof (replacing t+ by t- and reordering bounds in the integrals) works for 
the other mentioned resolution of the identity. According to Holder formula, we 
have: For any t & tt, 

\\(E t+ -E t )V(E t+ -E t )f\\l> = J* + V(fl {t ,r])(s) 2 ds 

<{r-tr+y 2 \\V{fl (t , r] )\\ cW -^ 
<c(r -^ +1 /2||/|| £2 . 

Thus for any e > 0, there exists rj > such that |7r| < n implies || (E t+ —E t )V(E t+ — 
Et)f\\c 2 < e - The proof is thus complete. □ 

The importance of strict causality lies in the next theorem we borrow from [11]. 

Theorem 3.4. The set of strictly causal operators coincides with the set of quasi- 
nilpotent operators, i.e., trace-class operators such that tracc(F™) = for any 
integer n > 1 . 

Moreover, we have the following stability theorem. 

Theorem 3.5. The set of strictly causal operators is a two-sided ideal in the set 
of causal operators. 

Definition 3.6. Let E be a resolution of the identity in the Hilbert space £ 2 ([0, T}; R"). 
Consider the filtration T E defined as 

Tf = <r{5 w (E x h), \<t,he C 2 }. 

An C 2 -valued random variable u is said to be T E -adapted if for any h e C 2 , the 
real valued process < E\u, h > is T E -adapted. We denote by E k (S)) the set of 
T E -adapted random variable belonging to Bp^f)). 

If E = (e\T, A e [0, 1]), the notion of T E adaptedness coincides with the usual 
one for the Brownian filtration and it is well known that a process u is adapted if 
and only if V^ v u(s) = for r > s. This result can be generalized to any resolution 
of the identity. 

Theorem 3.6 (Proposition 3.1 of [40]). Let u belongs to L P; i. Then u is T E - 
adapted if and only ifV w u is E-causal. 

We then have the following key theorem: 
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Theorem 3.7. Assume the resolution of the identity to be E = (e\T, A G [0, 1]) 
either E = (Id — e(i_^)r, A G [0, 1]). Assume that V is continuous from C 2 into CP 
for some p > 2 and that V is E-strictly causal. Let u be an element o/D| ; 1 (£ 2 ). 
Then, V\7 w u is of trace class and we have trace(V r V lv u) = 0. 

Proof. Since u is adapted, V w u is S-causal. According to Theorem 3.5, yV w ti is 
strictly causal and the result follows by Theorem 3.4. □ 

In what follows, E° is the resolution of the identity in the Hilbert space C? 
defined by e\f = /1[o,at] an d E° is the resolution of the identity denned by e\f = 
/1[(i-A)t,t]- The filtration J 7 ^' and are defined accordingly. Next lemma is 
immediate when V is given in the form V fit) — jl V(t, s)f(s) ds. Unfortunately 
such a representation as an integral operator is not always available. We give here 
an algebraic proof to emphasize the importance of causality. 

Lemma 3.8. Let V a map from £ 2 ([0, T]; R") into itself such thatV is E° '-causal. 
Then, the map T T Vfr T is E-causal. 

Proof. This is a purely algebraic lemma once we have noticed that 

(6) r T e r = (Id -er- r )T T for any < r < T. 
For, it suffices to write 

(7) T T e r f(s) = f(T - s)l[ , r ](T - s) 

= f(T - s)l [T -r,T\(s) = (Id-er-r)T T /(*), for any < s < T. 
We have to show that 

e r T T V^r T e r — e r T T V^T T or equivalently e r T T Vr T e r = T T Vr T e ri 
since e* — e r and r* = t t . Now, (7) yields 

Use (7) again to obtain 

eT~ r VT T e r — eT- r V(ld —eT- r )T T — {er- r V — eT- r VeT-r)T T = 0, 
since V is _E-causal. □ 

4. Stochastic integration with respect to Volterra processes 

In what follows, r\ belongs to (0, 1] and V is a linear operator. For any p > 2, 
we set 

Hypothesis I (p, 77). The linear map V is continuous from C p ([0, T]; R") into the 
Banach space Hol(ry). 

Definition 4.1. Assume that Hypothesis L(p, r\) holds. The Volterra process asso- 
ciated to V, denoted by W v is defined by 

W v {t) = S w {V(l [Qit] )), for all t G [0, T]. 



L. Dccrcuscfond 



9 



For any subdivision tt of [0, T], i.e., ir = {0 = to < t\ < . . . < t n = T}, of mesh 
7r|, we consider the Stratonovitch sums: 

. i rU+iAt 

(8) R 7r (t,u) = 6 Vf {J2g-l V-«(r)drl [tiit4+l) ) 



E i II V ^ru)(s) ds dr. 



* iG7r [t.At.t.+ lAt] 2 

Definition 4.2. We say that u is V- Stratonovitch integrable on [0,t] whenever the 
family R n (t,u), defined in (8), converges in probability as \ir\ goes to 0. In this case 
the limit will be denoted by f*u(s) o dW v (s). 

Example 1. The first example is the so-called Levy fractional Brownian motion of 
Hurst index H > 1/2, defined as 

rWTT75)i' (t - s)H " 1/2iB ' = 5(/ - I/2(1| °'' l,) - 

This amounts to say that V = I^I 1 ^ 2 . Thus Hypothesis I(p,H — 1/2 — 1/p) holds 
provided p(H — 1/2) > 1. 

Example 2. The other classical example is the fractional Brownian motion with 
stationary increments of Hurst index H > 1/2, which can be written as 

f K H (t,s) dB(s), 
Jo 

where 

r) = { -^^n\ -H,H-\,H + \,l- I )l m (r). 

The Gauss hyper-geometric function F(a, /?, 7, z) (see [24] J is the analytic contin- 
uation on C x C x C\{— 1, —2, . . .} x {z e C, Arg\l — z\ < ir} of the power series 



+00 



z , 



and 



T(a + k) 

(a) = 1 and (a) k = — = a(a + 1) . . . (a + k - 1). 

1 (a) 



We know from [33] that Kh is an isomorphism from £ p ([0, 1]) onto Z H+1 / 2 p an d 

K H f = I^ + x H -^C- 1/2 x^- H f. 
Consider Kh = I^ 1 o Kh- Then it is clear that 

f K H (t, s) dB(s) = I {K H )* T {l m ){s) dB(s), 
Jo Jo 

hence that we are in the framework of Definition ^.2 provided that we take V — 
(ICh)t- Hypothesis I(p, H — 1/2 — 1/p) is satisfied provided that p(H — 1/2) > 1. 
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Example 3. Beyond these two well known cases, we can investigate the case of 
multi-fractional Brownian motion defined as 

f K H(t) (t,s) dB{s), 
Jo 

for a deterministic function H. Consider the linear map Mh defined by 
M H : C 2 — ► C? 

/— > / K H{t) (t, s)f(s) ds. 
Jo 

The next result is extracted from [7] . 

Theorem 4.1. Let H belong to Hol(ry) for some r/ > 1/2 and be such that inf t H{t) > 
■q. Then Mh is continuous from C p into I* +1 / 2 p- 

If we consider Mr = Iq+ ° Mh, it holds that 

f K H{t) (t,s) dB{s) = f (Mff)5.(l [0 ,t])(«) dB(s). 
Jo Jo 

For, we remark that 

(M h )t = (Mh)t ° I T - = -(Mh)t o /q; 1 . 
Extending all the operators to distributions, we can write 

(-M*)r(l[o,t])(*) - {M H )* T e t {s) - f K H(r) (r,s) de t (r) = K H{t) (t,s). 

Jo 

Thus, Hypothesis I(p,r] — 1/2 — 1/p) is satisfied with V — (Mh)*t provided that 
p(»? - 1/2) > 1. 

Note that in this last example we do not have an expression of V as a kernel 
operator. This justifies the operator-theoretic approach of the sequel. 

The next theorem then follows from [7]. For the sake of completeness, we repro- 
duce its proof. 

Theorem 4.2. Assume that Hypothesis I(p, n) holds. Assume that u belongs to 
L p ,i- Then u is V -Stratonovitch integrable, there exists a process which we denote 
by D w Wu such that D w u belongs to L P (P ® ds) and 

(9) [ T U (s)o dW v {s)=5 w {Vu)+ f D w u{s) ds. 

Jo Jo 

Moreover, for any r <T, e r u is V Stratonovitch integrable and 

/ u(*)o dW v (s) = / (er«)(s)o dW v {s) =5 w (Ve r u)+ / D w u(s) ds. 
Jo Jo Jo 

Proof. Since u belongs to L Pi i, dP ® dr-a.s., the map (s i-> V^u s ) belongs to C p . 
Then, hypothesis l(p) entails that (s n> V(V™u) s ) is 77-Holdcr continuous. The 
map 

x [0, T] — > Hol(r/)* x Holfa) 

(w, r) — ^ ( £r , (« ^ V^(V» S )) 
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is measurable, hence the process 

D W U(C0, r) =< S r , V(V»(.) >Hol(r / )*,Hol W = V(V»(r), 

is measurable. Moreover, 



E 



f T \D w u(r)\ p 
Jo 



dr 



= E 



< E 



f I < e r , V(V»(.) >Hoi (ll )..Hoiw \ p dr 
Jo 



(10) 

Then, we have 



< c7^||u||£ p i . 



(11) E 



y- 



(V(V»(s) -£ w u(r)) d, (1/ 



[tiAt,t i+ iAt] 2 



< E 



< E 



< cE 



Z^ 0, 



(e a - £ r , T / (V^ V u)) H ol( J) )*,Hol(r ) ) 



ds dr 



[tiAt,*i + iAt] 2 



E^T // H e - - ^IIho1(,).II^(V»IILi(,) ^ dr 



[t;At,* i+ iAt] 2 



V-j- // (s-rf / |V^ M (a)| p dads dr 



< cItt^E 
<c|7rr|lu||? 



[tiAtjt.+ iAt] 

•t i+ iAt /-T 



E / / |V>a| p dadr 

t . en JUAt Jo 



Since 



E] J- D w u(r) ds dr = ^ £ w u(r) dr, 

*' £Tr ViAt,t i+ iAt] 2 

Eqn. (11) means that 

JJ V ^ru)(s) dr ds ^ w _ ^ C w U (r) dr. 



[tiAt,t i + 1 At] 2 



The remaining of the proof of (9) follows the classical proof for convergence of 
Stratonovitch sums as exposed in [28]. □ 

Theorem 4.3. Assume that Hypothesis I(p, r\) holds. Let u belong to L Pj i. If 
W w u is of trace class, then 

[ D w u(s) ds = tracc^V^it). 
Jo 
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Moreover, 

E [\trace(VV w u)\ p ] < c\\u\\l pi . 

Proof. Foreachfc, let (4>k,m, to = 1,-- • , 2 fc ) be the functions 4>k. m — 2 fc ^ 2 l[(„ l _i)2-'=.m2- fc )- 
Let Pfc be the projection onto the span of the <fik,m, since V w Vu is of trace class, 
we have (see [34]) 



trace(FV w p t w) = lim trace(P fc V V w p t u P k ). 

k— >+oo 



Now, 



tracc(P fe V\7 w uP k ) = (W w p t u, <j> k , m ® 4>k, m )c 2 ®c 2 

m— 1 

= Y, 2 l / V(V™u)(s) ds dr. 

m=l J (m-l)2- k At J (m-l)2- k At 

According to the proof of Theorem 4.2, the first part of the theorem follows. The 
second part is then a rewriting of (10). □ 

There is another result from [7] which is worth quoting for the sequel. 
Theorem 4.4. Assume that Hypothesis I(p, r\) holds. Let u belong to L Pj i. Then, 

| fu(s) o dW v (s)\f HoKn) 



(12) E 



where c does not depend on u. 



We can then follow the approach given for Stratonovitch integral as in [28] and 
show that we have a substitution formula. For p > 1, let T p be the set of random 
fields: 

u : R m — > L Pi i 

x i — > ((w, s) M> u(co, s, x)) 

equipped with the semi- norms, 

p K {u) = sup |Kar)|| L 1 
xeK 

for any compact K of R m . 

Corollary 4.5. Assume that Hypothesis I(p, rj) holds. Let {u(x), ,x E R" 1 } belong 
to T p . Let F be a random variable such that ((w, s) n> u(co, s, F)) belongs to L P; i. 

(13) [ u(s,F)o dW v (s)= [ u(s,x)o dW, 

Jo Jo 

Proof. Simple random fields of the form 

K 

u(u>, s, x) = Hi(x)ui(lj, s) 
i=i 



x=F 
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with Hi smooth and ui in L Pi i are dense in T p . In view of (12), it is sufficient to 
prove the result for such random fields. By linearity, we can reduce the proof to 
random fields of the form H(x)u(uj, s). Now for any partition ir, 

^ W (E J. \ H(F)V(u(u,.))(r) drl [ti , ti+l) ) 

= H(F)5™(Y1 T r +lA V(u(c,.))(r) dr l [t<it(+l) ) 



/■ti + lAt rti+iAt 

Y. / H'(F)V™FVu{r) ds dr. 

t . en JtiAt JtiAt 



On the other hand, 

V?(H(F)u(uj, r)) = H'(F)V7Fu(r), 

hence 



E^T // nV^(F)«)(«) d« dr 



** G7r [tiA*,t s+ iAt] 2 



** e,r [UAt,U+lAt] 2 



According to Theorem 4.2, Eqn. (13) is satisfied for simple random fields. □ 
Definition 4.3. For any < r < t < T, for u in h p , i, we define J r u(s) o dW v (s) 



as 



f u(s)o dW v (s) = f u{s)o dW v {s)- f u{s) o dW v {s) 

Jr Jo Jo 

= f e t u{s) dW v (s) - f e r u(s) o dW v (s) 
Jo Jo 

= S w {V{e t -e r )u) + f D w u{s) ds. 



Lemma 4.6. Let A and B be two continuous maps from £ 2 ([0, T]; R n ) into itself. 
Then, the map t t A <g> B (resp. At t <g> B) is of trace class if and only if the map 
A®t t B (resp. A®Bt t ) is of trace class. Moreover, in such a situation, 

trace(r T A <8> B) = trace(A t t B), resp. trace(Ar T <8> B) = trace(A <g> Bt t ). 

Proof. The map t t A ® B is of trace class if and only if for (h n , n > 1) a CONB of 

^2 | (t t A <g> B, h n ®h n )cA < 00 . 
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But, 

\ (t t A <8> B, h n ® h n )cA = ^ \(r T Ah n , Bh n )&\ 

n n 

= \( Afl n, r T Bh n ) C 2\ 

n 

= y^|(A® t t B, h n <g> /i„)£2|, 

n 

since t t is self-adjoint in £ 2 . The first result follows. The second part follows by 
adjunction. □ 

Corollary 4.7. Let u e L 2i i such that \7 W ® r T T/u awd V 1 *"® Vt t u are o/ irace 
c/ass. T/ien, t t V w ® T^u and V w t t ® are of trace class. Moreover, we have: 

trace(V vv Or T yu) = trace(r T V lv <g) Vu) 

and trace(V lv (8) (Vt t )u) = trace(V M/ r T (g) 
Proof. For u simple, i.e., of the form 

n 

u (u, s) = ^2Ui(uj)gi(s), 

3 = 1 

the result follows from Lemma 4.6. Such random fields are dense in F p and according 
to Theorem 4.3, the trace function is continuous on T p hence the result is satisfied 
in full generality. □ 

Theorem 4.8. Assume that Hypothesis I(p, n) holds. Let u belong to L Pi i and let 
Vr = T t VTT- Assume furthermore that V is E - causal and that u = u o Q^ 1 is 
T E ° -adapted. Then, 

r T-r ft 

(14) / t t u{s)o dW v (s) = / V T (l [r .t]u)(s) dB T (s), < r < t < T, 

JT~t Jr 

where the last integral is an ltd integral with respect to the time reversed Brownian 
motion B T (s) = B(T) - B(T - s) = Q T (B)(s). 

Proof. We first study the divergence term. In view of 3.2, we have 

8 B (V(eT- r — e-T-t)T T u o 9 T ) = S B (VT T (e t — e r )u o T ) 

= S B (r T Vr(e t — e r )u o © T ) 

= 5(V T (e t - e r )u)(&) 



= J V T (l [r , t] u)(s) dB T (s). 



According to Theorem 3.8, (Vt)* is E causal and according to 3.3, it is strictly 
Eo causal. Thus, Theorem 3.7 implies that VV(et — e r )u is of trace class and 
quasi-nilpotent. Hence Lemma 4.7 induces that 

t t Vtt t ® T T Vr T (et — e r )u 

is trace-class and quasi-nilpotent. Now, according to Theorem 3.1, we have 

t t Vtt t <g> T T Vr T (e f — e r )u = V(Vr T (eT- r — er-t)u o @t). 

According to Theorem 4.2, we have proved (14). □ 
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Remark 4.1. Note that at a formal level, we could have an easy proof of this 
theorem. For instance, consider the Levy fBm, a simple computations shows that 
Vt = Iq-^ 1 ^ 2 I or an U T. Thus, we are led to compute trace(/ ( ^~ 1/ ' 2 VM). If we had 
sufficient regularity, we could write 

trace(/£~ 1/2 Vu) = f [ (s - r)^ 3 / 2 V s u(r) dr ds = 0, 
Jo io 

since V s u(r) = for s > r for u adapted. Obviously, there are many flaws in these 
lines of proof : The operator 1 ' 2 V« is not regular enough for such an expression 
of the trace to be true. Even more, there is absolutely no reason for VtVu to be 
a kernel operator so we can't hope such a formula. These are the reasons that we 
need to work with operators and not with kernels. 

5. VOLTERRA DRIVEN SDES 

Let (5 the group of homeomorphisms of R™ equipped with the distance: We 
introduce a distance d on © by 

d(<P, 4>) = 4>) + pi'P' 1 , 

where 

\<p(x) - (j>{x)\ 



l+SUp M< , 



, .•supi^i^l^a;)-^)! 
Then, © is a complete topological group. Consider the equations 

(A) X r>t =x+ I a{X riS ) o dW v {s), < r < t < T. 

J r 

(B) Y r ,t=x-J a{Y Sit )o dW v (s), < r < t < T. 

Definition 5.1. By a solution of (A), we mean a measurable map 
n x [0,T] x [0,T] — > © 

(ui, r, t) i — > (x i ^ X rit (u,x)) 
such that the following properties are satisfied : 

(1) For anyQ <r <t < T, for any x G R", X rtt (u, x) is a{W v (s), r < s < t}- 
measurable, 

(2) For any < r < T , for any x G R", the processes (u, t) i— > X r t {uj, x) and 
(cj, t) i ^ X~ t (u), x) belong to L Pj i for some p > 2. 

(3) For any < r < s < t, for any x G R", the following identity is satisfied: 

X rtt (uj, x) = X Syt (uj, X r ^ s (uj, x)). 

(4) Equation (A) is satisfied for any < r < t < T P-a.s.. 
Definition 5.2. By a solution of (B), we mean a measurable map 

ft x [0, T] x [0, T] — > © 

(lj, r, t) i — > (x n> Y rtt (u),x)) 
such that the following properties are satisfied : 

(1) Forany0<r<t< T, for any x G R n , Y r>t (u}, x) is a{W v (s), r < s < t}- 
measurable, 
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(2) For any < r < T , for any x € R™, the processes (u>, r) i-> Y rt t(u>, x) and 
(w, r) i-» Y~ t (u, x) belong to L Pj i for some p > 2. 

(3) Equation (B) is satisfied for any < r < t < T P-a.s.. 

(4) For any < r < s < t, for any x € R" ; the following identity is satisfied: 

Y rt t(u), x) = Y r , s (u>, F s ,t(w, x)). 
At last consider the equation, for any < r < t < T, 

(C) Z rtt = x- [ V T {ooZ.+ l M )(s) AB t {s) 

where B is a standard n-dimensional Brownian motion. 

Definition 5.3. By a solution of (C) ; we mean a measurable map 

fl x [0,T] x [0,T] — > © 

(w, r, t) i — > (x Z rit (w,a;)) 

swc/i i/iai t/ie following properties are satisfied : 

(1) For anyO<r <t< T, for any x £ R", Z r , t {u, x) is a{B T (s), s<r <t}- 
measurable, 

(2) For any < r < t < T, for any x e R", the processes (u, r) n> Z r _ t (cj, x) 
and (u>, r) H> Z~1(uj, x) belong to L p a for some p > 2. 

(3) Equation (C) is satisfied for any < r < t < T P-a.s.. 

Theorem 5.1. Assume that Vt is an E causal map continuous from CP into I a , p 
for a > and p > 4 such that ap > 1. Assume a is Lipschitz continuous and 
sub-linear, see Eqn. (21) for the definition. Then, there exists a unique solution to 
equation (C). Let Z denote this solution. For any (r, r'), 

E[\Z riT -Z r ,, T f] < c \r-r'r. 

Moreover, 

(co, r) i V Z r ^ s (uj, Z Sit (uj, x)) S L p ,i, for any r < s < t < T. 

Since this proof needs several lemmas, we defer it to Section 6. 

Theorem 5.2. Assume that Vt is an E° causal map continuous from C p into I a _ p 
for a > and p > 2 such that ap > 1. For fixed T, there exists a bisection between 
the space of solutions of Equation (B) on [0, T] and the set of solutions of Equation 
(C). 

Proof. Set 

Z r ^ T (uj, x) = y T _ I . iT (9y 1 (tI'), x) 

or equivalently 

(15) Y r ,T{u, x) = Z T _ r! T(e T (w), x). 

According to Theorem 4.8, Y is satisfies (B) if and only if Z satisfies (C). The 
regularity properties are immediate since £ p is stable by r T . □ 

The first part of the next result is then immediate. 
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Corollary 5.3. Assume that Vt is an E° causal map continuous from CP into la.p 
for a > and p > 2 such that ap > 1. Then Equation (B) has one and only solution 
and for any < r < s < t, for any x e R™, the following identity is satisfied: 

(16) Y rit (cj, x) = y r , s (w, y S)t (w, xj). 

Proof. According to Theorem 5.2 and 5.1, (B) has at most one solution since (C) 
has a unique solution. As to the existence, point (1) to (3) are immediatly deduced 
from the corresponding properties of Z and Equation (15). 

According to Theorem 5.1, (w, r) h-» Y r ^ s (oj, Y s j(oj, x)) belongs to L Pi i hence we 
can apply the subsitution formula and we get: 



x=Y Sit (uj, x) 



(17) Y r , B (u>, y 8it (w, x)) = Y Stt (w, x) - f a(Y T , s (u, x)) o dW v (r) 

J r 

= x - f a{Y Tjt {uj, x) o dW v (r) 

•/ s 

- j' a(Y TtS (w, Y Stt (u, x))) o dW v (r). 



Rr.t — 



Set 

Y Tt t(u), x) for s < r < t 

y T ,«(w, y s ,t(w, a;)) for r < r < s. 

Then, in view of (17), i? appears to be the unique solution (B) and thus R Si t(oJ, x) — 
Y 8>t (tJ, x). Point (4) is thus proved. □ 

Corollary 5.4. For x fixed, the random field (Y r j(x), < r < t < T) admits a 
continuous version. Moreover, 

E [\Y r , s (x) - Y r ,^(x)\P] < c{\ + \x\P)(W - s\ pr > + \r - r'D. 

We still denote by Y this continuous version. 

Proof. W.l.o.g. assume that s < s' and remark that F S;S /( X ) thus belongs to 
a{B^, u > s}. 

E[\Y r , s (x)-Y r , iS ,(x)\ p ] 

< c(E [\Y riS (x) - Y r , iS (x)\ p ] + E [\Y r , iS (x) - Y r , , s , {x)\ p \) 
= c (E [\Y r ^(x) - Y r , iS (x)\ p ] + E [\Y r , iS (x) - Y r ,, s (Y s , s ,(x))\ p )) 
= c (E [\Z s _ riS (x) - Z s _ r , iS (x)\ p ] + E [\Z s _ r ,, s {x) - Z s . r ,, s (Y s , s ,(x))\ p }) . 
According to Theorem 6.2, 

(18) E [\Z s _ r , s (x) - Z s _ r ,, s (x)\ p ] < c\r ~ rT(l + \x\ p ). 

In view of Theorem 4.8, the stochastic integral which appears in Equation (C) is 
also a Stratonovitch integral hence we can apply the substitution formula and say 

Zs-r',s(Ys,s'(x)) = Z s _ r i ,s(y)\ y —Y a s ,( x ) ' 

Thus we can apply Theorem 6.2 and we obtain 

E[\Z a -r,, a (x) - Z s _ r ,, s {Y s , 8 ,(x))\ p ] < cE[\x -Y SiS ,(x)\ p ] . 
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The right hand side of this equation is in turn equal to E [\Zo, S ' — Z s i- SiS i(x)\ p ] , 
thus, we get 

(19) E [\Z s _ r ,, s (x) - Z s _ r ,, s (Y s , s ,(x))f] < c(l + \x\p)W - s r 
Combining (18) and (19) gives 

E [\Y r , a {x) - Y r ,, s ,(x)\ p ] < c(l + \x\ p )(\s' - + \r- r'D, 

hence the result. □ 

Thus, we have the main result of this paper. 

Theorem 5.5. Assume that Vt is an E causal map continuous from C p into I a ^ p 
for a > and p > 4 such that ap > 1. Then Equation (A) has one and only 
solution. 

Proof. Under the hypothesis, we know that Equation (B) has a unique solution 
which satisfies (16). By definition of a solution of (B), the process Y^ 1 : (w, s) i— > 
Y s 7 1 (w, x) belongs to L Pj i hence we can apply the substitution formula. Following 
the lines of proof of the previous theorem, we see that Y^ 1 is a solution of (A). 

In the reverse direction, two distinct solutions of (A) would give raise to two 
solutions of (B) by the same principles. Since this is definitely impossible in view 
of Theorem 5.3, Equation (A) has at most one solution. □ 

6. The forward equation 

Lemma 6.1. Assume that Hypothesis I holds and that a is Lipschitz continuous. 
Then, for any < a < b < T, the map 

V T oo- : C([0,T], R")— ►C([0,7 , ] ) Il n ) 

0i — >V T (aoip 1 [0)6] ) 

is Lipschitz continuous and Gateau differentiable. Its differential is given by: 

(20) dV T o a(<p)[ip] = V t (<j' °<t> # 
Assume furthermore that a is sub-linear, i.e., 

(21) \a(x)\ < c(l + |.t|), for any x G R". 
Then, for any V> £ C([0, T], R n ), for any t e [0, T], 

\V T {o- o ^)(t)\ < cT" +1 / p (l+ / \tp(s)\ p ds) 

Jo 

< Cl">+Vp(i + Woo ). 

Proof. Let tp and <fr be two continuous functions, since C([0, T], R") is continuously 
embedded in C p , Vt(o- o ip — a o <f>) belongs to Hol(ry). Morevover, 

sup \V T (a o V l [0)6] )(i) - Vt(t o l [a . fc] )(t)| < c ||Vt((ct o V - cr o 0) 1 [a , b] ) || H oi(r,) 

< c||(cto^ -ao<f>) l [a , fc] ||£P 

<c\\<p-1p\\ C P([a,b}) 

<csup|V(i)-0(i)|, 

since cr is Lipschitz continuous. 
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Let if; and ip two continuous functions on [0, T]. Since a is Lipschitz continuous, 
we have 

a(ip(t) + e<t>{t)) = a(ip(t)) + e [ a'(uip{t) + (1 - u)4>{t)) du. 



Moreover, since a is Lipschitz, a' is bounded and 



a'(uip(t) + (1 - u)4>(t)) du 



dt < cT. 



This means that (t i-> J cr / (uV'(t) + (l — u)(j>(t)) du) belongs to CP . Hence, according 
to Hypothesis I, 

\\V T ( [ <t\u^{.) + (1 - u)<j>{.)) du)||c < cT. 

Thus, 



o 



lim £ 1 (Vt(ct o M) + e4>)) — Vr(cr o tp j) exists, 

and Vt ° cr is Gateaux differentiablc and its differential is given by (20). 
Since cr o ifj belongs to C([0,T], R"), according to Hypothesis I, we have: 

/ f t \ Vp 

\Vt(<to ^){t)\ < c (j s r iP\a(^(s))\ p ds 

/ft \ Vp 

<cW / (l + |^(*)| p )ds 

< cT" +1 /f(l + I^IU). 

The proof is thus complete. □ 

Following [38], we then have the following non trivial result. 

Theorem 6.2. Assume that Hypothesis I holds and that a is Lipschitz continuous. 
Then, there exists one and only one measurable map from x [0, T] x [0, T] into <S 
which satisfies the first two points of Definition (C). Moreover, 

E [\Z r , t {x) - Z r , it (x')\ p ] < c(l + \x\ p V 1x7) (|r - r'| p " + |z - x'\ p ) . 

Note even if x and x' are replaced by a{B T (u), t < u} measurable random variables, 
the last estimate still holds. 

Proof. Existence, uniqueness and homeomorphy of a solution of (C) follow from 
[38] . The regularity with respect to r and x is obtained as usual by BDG inequality 
and Gronwall Lemma. For x or x' random, use the independance of a{B T (u), t < 
u} and a{B T {u), r A r' < u < t}. □ 

Corollary 6.3. Assume that Hypothesis I holds and that a is Lipschitz continuous 
and sub-linear. Let Z be a solution of (C). Then, for any x G R", for any 
< r < t < T, we have 



E[\Z rit (x)\"] <c(l + |af)e 
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Proof. According to BDG inequality and Lemma 6.1, we have: 



E[K t (z)H <c(|af + E 



f \V T {o-oZ., t {x)l [rt 

J r 



<cT w+1 (l + |a;| p + E 
< cT^ +1 (l + \x\p + jf E [\Z Tjt \P] dr) 



Z Tt (x)| p dr ds 



) 



We then conclude by Gronwall Lemma. 



□ 



Theorem 6.4. Assume that Hypothesis I holds and that a is Lipschitz continuous 
and sub-linear. Then, for any x G R", for any 0<r<s<t<T, (u>, r) <— > 
Z r _ s (u, Z s _t{x)) and (co, r) i-> Z~\{uj, x) belong to L Pi i. 

Proof. According to [18, Theorem 3.1], the differentiability of to M> Z rtt (oj, x) is 
ensured. Furthermore, 

V„Z r , t = -V T (aoZ. t tl lr ,t])(u) - I VT{v'(Z.,t)-VuZ., t l [r , t ])(s) dB(s), 

J r 

where a' is the differential of a. Since Vt is continuous from C p in to itself and a 
is Lipschitz, according to BDG inequality for r < u, 



E[|V„Z r , t H 



< cE 



|VHcT0Z., t l M ])(«)| P 



+ cE 



< c 1 + E 



|Z Tt | p dr dw 



/ \V T (v'(Z.,t)VuZ., t l[r,t])(s)\ p ds 

J r 

V u Z T j\ p dr ds 



E 



<c(l + F, \Z T . t \ v (t m+1 -T m+1 ) dr +E jf |V u ^,t| p (t pT?+1 -r p " +1 ) dr ^ 



< ct m+1 1 + E 



/' 



\Z T ,t\ p dr 



+ E 



f |v u z T>t |' 

J r 



dr 



Then, Gronwall Lemma entails that 

E[|V„Z r;t |f] < c ( 1 + E 



/Vr,t| P dr 

J r 



The integrability of E [|V u Z r! i| p ] with respect to u follows. 

Now, since 0<r<s<t<T, Z Sit (x) is independant of Z riS (x), thus the 
previous computations still hold and (lo, r) h-> Z r _ s {uj, Z s _ t (x)) belong to L Pj i. 

According to [35], to prove that Z~ t (x) belongs to B Pj i, we need to prove 

(1) for every h 6 C 2 , there exists an absolutely continuous version of the process 
{t i y Z~l{Lo + th, x)), 

(2) there exists DZ~l , an £ 2 -valued random variable such that for every h 6 
C 2 , 

- (Z-j (w + th, x) - Z~l (w, x)) ^ f DZ~l (s)h(s) ds, 
1 Jo 

where the convergence holds in probability, 
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(3) DZ~l belongs to C 2 (Q, C 2 ). 
Wc first show that 



(22) 
Since 



E 



dZ r1 



dx 



(w, Z-^x)) 



is finite. 



dZ, 



r, t 



(w, x) = Id + jf* Vr(a'(Z., t (x))^^)( S ) dR(s), 



OT J r OX 

Ito's formula and BDG inequality yield, for any q e R, for any 0<r<u<t<T, 



E 



sup\d x Z Utt (x)\ 2 " 

U<V 

< c + cE 
+ cE 



Mhr,t]cr\Z.4x))d x Z.4x))(s) 2 \d x Z s , t (x)\ 2 ^ ds 
\d x Z s ^x)\^ 2 V T {l [rA a\Z.4x))d x Z.4x)){s) 2 ds 



Let & v = swp u<v \d x Z Utt (x)\, since Vt maps £ p to we have 
rt / r s \ 2/p 



E[Q 2 v q ] <c + cE 



/t / rs \ -VP 



+ cE 



er 2 



|a a Z T>t (a;)n dr 



2/pN 



ds 



Hence, 



E[e^]< c (i + £E[e 2 «] ds), 



and (22) follows by Gronwall Lemma. Since Z r> t{ui, Z~ t (cj, x)) = x, the implicit 
function theorem imply that Z~\(x) satisfies the first two properties and that 
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VZ r ^ t (Lo, Z~J(x)) + -^(u>, Z-j(x))VZ-j(uj, x). 
It follows by Holder inequality and Equation (22) that 

\\DZ-j(x))\\ pA < c||Z I . it (x))||2 Pi i||(a a! Z I . it (a:))- 1 ||2 P) 
hence Z~\ belongs to L Pi i. 
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